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In a covariant model where constituent quarks and diquarks interact through quark 
exchange, the Bethe-Salpeter equation in ladder approximation for octet and de- 
cuplet baryons is solved. Quark and diquark confinement is thereby effectively 
parametrised by choosing appropriately modified propagators. Numerical results 
for the baryon masses are presented. In a second step electromagnetic, weak and pi- 
onic currents are coupled to the bound state according to Mandelstam's technique. 
The arising matrix elements are evaluated in a generalised impulse approximation 
and observable form factors are extracted. 

1 Introduction 

With the advent of a new generation of continuous beam facihties Hke TJNAF, 
COSY and ELSA which investigate baryon observables at a momentum regime 
of = 1 — 5 GeV to a very high precision, a new era of hadron phenomenology 
starts. It is well known, that neither pure hadronic models nor non-relativistic 
quark models are able to explain the wealth of phenomena arising in the transi- 
tion region between non-perturbative strong interaction physics and the regime 
where QCD perturbation theory sets in. Thus, it proves helpful to look for 
effective but also efficient degrees of freedom within a covariant approach. 

Inspired by studies in the hadronized Nambu Jpna-Lasinio modelld and the 
Global Color model □, in two recent publications we worked out a covariant 
description of baryons, where quarks and diquarks enter as effective degrees of 
freedom. In the past, the notion of diquarks, taking into account some of the 
unknown features of non-perturbative QCD, served as a valuable tool for vari- 
ous hadron phenomenal. Apart from being manifestly covariant, our approach 
includes an effective parametrisation of quark and diquark confinement. The 
framework to examine relativistic bound states is given by describing the quark 
exchange between quarks and diquarks within the Bethe-Salpeter approach in 
ladder approximation. In this way bound state masses and wave functions are 
determined. 
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2 Bethe-Salpeter Equation for Baryons 



In the covariant diquark-quark modelcla, a baryon is described as a bound state 
of a constituent quark and a scalar (0^) or axialvector (1^) diquark interacting 
through quark exchange. The corresponding Bethe-Salpeter equation (BSE) 
for octet and decuplet baryons in ladder approximation is given by 0, 

Wj^[-V^§j5Si-q)r -^-^r'si-q)r ) Wip';p) 

^•^Sip-^P) = -2\gXSiPa)D-^nPb)J0^l'S{^qh^-fiSip';P)- (2) 

This equation determines the matrix-valued Bethe-Salpeter wave functions of 
octet baryons 'i'sip', P) and 5'g(p; P) and decuplet baryons ^'^q (p; P), respec- 
tively which are to be projected on positive parity and spin 1/2 or spin 3/2. In 
flavor space the wave functions correspond to pure octet and decuplet states. 
They explicitly depend on the total momentum P of the bound state and the 
relative momentum p' orp between the two constituents. In eqs.(|l|,^ the quark 
propagator is denoted by S{q), the scalar diquark propagator by D{q) and the 
axialvector diquark propagator by D^'^(q), see below. 



2. 1 Parametrisation of confined constituents 

One of the striking features of QCD is the absence of particles carrying a net 
color charge in the physical spectrum. Accordingly, the decay of a baryon 
into free quarks, or in our approach quarks and diquarks, is forbidden. Note, 
in a diquark-quark approach, diquarks as constituents of color-less baryons 
necessarily form an anti-triplet representation of the color group. 

From Dvson-Schwinger studies of QCD, assuming a suitable quark-quark 
interaction eI, it is known how quark confinement may be realized: The quark 
self-energies appearing in the propagator 

S{p) = -ij.pa,{p) + as{p) = [n-pA{p) + (3) 

have no poles on the real p^— axis. This feature can be parametrised by the 

"We use an Euclidean space formulation with {"/ii,"/!^} = 2<5^i/,7p = "/i^ and 
P-9 = X]Ll =Pi1i- 
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choice 

l-expM(l + ^) 

cr„(p) = — — and as{p) = mqa.o{p). (4) 

+ 

Obviously the exponentials remove the mass poles, however, thereby introduc- 
ing an essential singularity at = — oo. 

In refsLrQ it has been shown that confined diquarks can be obtained when 
one uses an appropriate irreducible quark-quark kernel (beyond ladder approx- 
imation) in the diquark Bethe-Salpeter equation. A useful parametrisation of 
confined diquarks is now given by the choice 

Div) - D^^ip) = (5) 

p'^ + rrig p^ + 



where 



p 

= I - cxp(-(i(H ^)), mrf = TTis or ma- (6) 

Note, for point-like (tree level) quarks and diquarks the exponentials are ab- 
sent. The parameter d basically controls the damping of the confining propaga- 
tors for space-like momenta. As shown below, this prescription of confinement 
excludes the decay of a baryon into a free quark and diquark. 
In order to obtain a finite quark-diquark interaction in momentum space and 
to simulate the intrinsic structure of diquarks in a crude way, we modify the 
propagator of the exchanged quark in the BSE (|^,§) according to (see also 
reffl) 

S{q) ^ ~S{q) = S{q) (^J^) ' 
2.2 Numerical solution and results 
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Before solving the BSE (|l|,y) we project the wave functions ^'8(p; i-'),5'g(p; P), 
'i'iQ{p;P) onto the relevant spin and parity content, restrict ourselves to so- 
lutions with positive energjLand decompose the resulting covariant Lorentz 
tensors in the Dirac algebrau'u. Furthermore, the scalar functions multijobying 
each tensor are expanded in Chebyshev polynomials of the second kindlj. In 
this way, the BSE is transformed into a system of coupled hijmogeneous integral 
equations, which can be solved by an iterative algorithniEj. The eigenvalue of 
the integral equation corresponds to the diquark-quark coupling strength as a 
function of the bound state mass M; the eigenvectors are used to reconstruct 
the baryon wave functions. In the actual calculation SU(3)flavor is broken by 
the strange quark mass mst and isospin symmetry is assumed. 
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0.5 
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rUst (GeV) 




0.66 


0.63 






1 


0.73 


9a 




10.35 


10.89 


9s 




11.99 


10.77 


Ma (GeV) 


1.116 


1.125 


1.120 


Me (GeV) 


1.193 


1.163 


1.154 


Mh (GeV) 


1.315 


1.325 


1.314 


Me. (GeV) 


1.384 


1.383 


1.373 


Mh. (GeV) 


1.530 


1.566 


1.549 


Mo (GeV) 


1.672 


1.723 


1.699 






0.0022 


0.0016 



Table 1: Octet and decuplet masses for two different parameter sets (I), (11). 

In table (|l|) we display our results for the octet and decuplet masses ob- 
tained for two different parameter sets □. The nucleon and delta mass have 
been used as input to determine the coupling constants ga and gs and are 
therefore not displayed. Here, ^ is defined by the ratio ^ = mab/ima + nT'b), 
where ab is the flavor content of a diquark. As can be seen, it is possible to get 
an excellent overall fit to the experimental spectrum, as indicated by the given 
X^. Note that the octet-decuplet mass differences are solely due to the differ- 
ent interaction in the octet and decuplet channels; the scalar and axialvector 
diquark masses are assumed to be equal. The main advantage of our approach 
is, however, that we furthermore have access to the relativistic Bethe-Salpeter 
wave functions, which can be used in calculations of form factors and other 
non-static observables, see below. 

The effect of using confining propagators is clearly seen in figure (Q). 
When using tree level (non-confining) propagators we are able to determine 
the eigenvalue of the integral equation up to the diquark-quark threshold, here 
at M/rUq — 2. At threshold, the bound state decays into a free quark and 
diquark, which is, of course, unphysical. In case of confining propagators, the 
threshold vanishes and the baryon masses are allowed to be larger than the 
sum of quark and diquark masses. 
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Figure 1: Eigenvalue of the nucleon BSE versus bound state mass M . Threshold effects are 
absent for confining propagators. 



3 Nucleon form factors 

In the following we report on a calculation of nucleon form factors within the 
covariant and confining diquark-quark model where only scalar diquarks with 
nis — rriq are considered El. 

In the Bethe-Salpeter approactLform factors are most conveniently calcu- 
lated via Mandelstam's formalismllil: 

Besides the normalised BS wave function 5*, the 5-point function Tq enters, de- 
scribing the coupling of the probing current to all internal lines of the diquark- 
quark system. 

In the following we use a generalised impulse approximation, where only 
the coupling to the quark and to the scalar diquark is considered, 

T^ip',P';p,P)=Tl{p',P';p,P)+T'i{p',P';p,P), (9) 

thereby neglecting the coupling to the exchanged quark. Using this approxi- 
mation, the matrix elements split into a quark and a diquark part, given by 
the loop integrals, containing the fully amputated BS vertex function x{P',p), 

^(^^^^'^ ^^mI ^m^Pf)S{p+)Tl)iP+,P^)Sip^)Dip,)x{Pr,P^) 

(10) 
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and 

□ (11) 

For details see rem. 

3. 1 Electromagnetic form factors 

To calculate electromagnetic (e.m.) form factors, the quark-photon and the 
diquark-photon vertex functions entering (^0|) and have to be specified. 
However, their choice can be restricted using gauge symmetry of QED, mani- 
fest in the Ward and Ward-Takahashi identity. 

A longitudinal quark-photmi vertex function in agreement with gauge symme- 
try is the Ball-Chiu vertexllj, 



-z(i?(p^)-B(fc^))}, (12) 

which is fully determined by the self-energy functions A{p'^) and B{p^) of the 
quark propagator (^. 

As has been shown in reftj a longitudinal vertex function of a scalar di- 
quark which guarantees gauge invariance at the constituent level is given by 



Tiip, k) = -ip + fc) /"^^Lp"'^^^ ' ^13) 



which is fully determined by the self-energy function of the diquark propagator. 

In the generalised impulse approximation, the electromagnetic nucleon cur- 
rent is denoted by 

MQ^) = 0, J^(Q') + QdJ^,iQ% (14) 

where the quark and diquark charges, Qq and Qd, appear. Using the familiar 
decomposition of the on-shell e.m. nucleon current, where A+ is the projector 
onto positive energy states, 

J^(g^) = A+{Pf)[iMB{F,{Q^) - F„,{Q'))^ + F,n{Q'h^]A+{P^). (15) 

the electric and magnetic form factors, Fe{Q^) and FmiQ'^), can be extracted 
in a straightforward way. The left hand side of this equation is calculated with 
the loop integrals diq,nih. 



3.2 Pion-nucleon and weak axial form factor 



Considering only scalar diquarks, the coupling of an external pion or an axi- 
alvector current to the nucleon is especially simple: Due to parity these currents 
solely interact with the quark but not with the scalar diquark. Furthermore, 
the quark-pion vertex (= pion Bethe-Salpeter amplitude) in the chiral limit 
(m^r = 0) is given by the scalar quark self-energy 

^g(p2^0,/)=75^r^ (16) 

as an immediate consequence of chiral symmetry. The pion-nucleon matrix 
element is commonly parametrised by 

JI{Q^) = K+{Pf)[^,g.NN{Q^K]^+{P^). (17) 

including the pion-nucleon form factor g-r^NNiQ'^)- The left hand side of this 
equation is calculated with eq. (^) supplemented with the quark-pion vertex 
&■ 

The coupling of an axialvector current to the quark is implemented via the 
vertex function 

^5f.iP^k) = [i 7^ (18) 



p2 _ 



{A{p^) - A(fc2))2Z±^ - ^{B{p^) + B{e)) 



75y, 



which was constructed to obey the chiral Ward identitjQ. The invariant func- 
tions entering this expression are again determined by the self-energy functions 
of the quark propagator. Using the decomposition 

j;{Q^) - A+(P;)y [7^5^(Q2) + Q^g^(Q2)]^^A+(P0 (19) 

of the nucleon axialvector current, the form factor gA{Q'^) can be extracted. 



3.3 Results 

In table (||)j-.we present our results for the static form factors {Q^ — 0) of 
the nucleonl3 obtained with tree level and confining propagators. Due to the 
employed generalised impulse approximation, charge conservation is violated 
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Figure 2: Electromagnetic nucleon form factors. The nucleon mass in all cases is chosen 
to be M = 1.9mg and the amplitudes entering the calculation of the matrix elements are 

obtained with A = 2mq , 

by a few percent Furthermore, the magnetic moments of proton and neutron 
are too smah, adiich was also observed in similar approaches including only 
scalar diquarksll3'tZH13. The pion- nucleon form factor and gA, however, are in 
a reasonable agreement with the experimental values. This indicates the strong 
coupling of the corresponding currents to the quark due to chiral symmetry. 

The momentum dependence of the e.m. form factors can be seen in fig- 
ure (||). In the calculation with confining propagators and d — 10, we get a 
Q^-behaviour of the form factors similar to the calculation with tree level prop- 

''Note, ther4riclusion of the diagram, where the photon also couples to the exchanged 
quark leads tdij Qp + Qn = 1- 
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Nucleon form factors 




Tree level prop. 


Conf. prop. , d = 10 


Exp. 


Qp 


1 + 4.1 • 10^3 


1 + 4.7-10-3 


1 


Qn 


-8.2- 10-3 


-9.4-10-3 





Hp 


1.32 


1.34 


2.79 




-0.64 


-0.65 


-1.91 


9A 


1.39 


1.41 


1.25 


9nNN 


10.59 


10.89 


13-14.5 



Table 2: Numerical results for the static nucleon form factors. The magnetic moments are 

given in units ofe/2M. 

agators. However, when using confining propagators we observe that the slopes 
arc less sensitive to the model parameters and to the effects of the diquark- 
quark pseudo-threshold. In the figure showing the electric proton form factor 
(-F(f (Q^)) we also display the empirical dipole fit. It is clearly seen that our 
charge radii are smaller than the experimental value. However, in the present 
stage of the investigation it makes no sense to tune the various parameters 
of the model to get closer to the dipole curve and to reproduce the charge 
radius of the proton. This has to be done, after axialvector diquarks and the 
coupling to the exchanged quark are included in the form factor calculation. 
But even without fine-tuning we obtain qualitatively acceptable results which 
indicates that the gross features of baryon structure can be described within 
our approach. 

4 Conclusion 

In this talk we have presented a covariant diquark-quark bound state approach 
for baryons where the constituents interact through quark exchange. Using an 
effective parametrisation of confinement we calculated with a Bethe-Salpeter 
equation in ladder approximation masses and wave functions of octet and de- 
cuplet baryons. An excellent description of the baryon spectrum is obtained 
with the assumption that SU(3)flavor is broken by the strange quark mass. 

As a first application of this baryon model, we investigated ele(;tromag- 
netic, pionic and axial form factors of the nucleon, considering only scalar 
diquarks. While for a complete and quantitatively reasonable description of 
these observables axialvector diquarks are obviously necessary, we neverthe- 
less conclude from our studies, that this approach provides a sound basis for 
baryon phenomenology. 
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The advantages of the covariant and confining diquark-quark model will 

be even more visible, when applied to medium energy reactions. Currently, 
electromagnetic and hadronic A-production and virtual Compton scattering 
are therefore under investigation. 
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